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Abstract 

We present a new perturbative expansion for pionless effective field theory with Coulomb interac¬ 
tions in which at leading order the spin-singlet nucleon-nucleon channels are taken in the unitarity 
limit. Presenting results up to next-to-leading order for the Phillips line and the neutron-deuteron 
doublet-channel phase shift, we find that a perturbative expansion in the inverse 1 S'o scattering 
lengths converges rapidly. Using a new systematic treatment of the proton-proton sector that 
isolates the divergence due to one-photon exchange, we renormalize the corresponding contribu¬ 
tion to the 3 H- 3 He binding energy splitting and demonstrate that the Coulomb force in pionless 
EFT is a completely perturbative effect in the trinucleon bound-state regime. In our new expan¬ 
sion, the leading order is exactly isospin-symmetric. At next-to-leading order, we include isospin 
breaking via the Coulomb force and two-body scattering lengths, and find for the energy splitting 
(£ B ( 3 He) - E b ( 3 H)) nlo = (-0.86 ± 0.17) MeV. 
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I. INTRODUCTION 


It has been known for a long time that two-nucleon (NN) scattering at very low energies 
can be described by the effective range expansion [T], from which deuteron and 1 S'o virtual 
state properties emerge without information about details of the strong interaction [2], At 
these energies effects from pion-exchange physics cannot be resolved. It is therefore possible 
to describe such a system using only nonrelativistic nucleons as degrees of freedom that 
interact via short-range (contact) forces [3bT2]. The systematic approach to implement this 
procedure, known as pionless effective field theory (pionless EFT), is based on the experi¬ 
mental fact that the NN S -wave scattering lengths are much larger than the corresponding 
effective ranges, so that a nonperturbative resummation of non-derivative two-body con¬ 
tact interactions is required at leading order (LO) to reproduce the shallow NN bound and 
virtual states. 

The extension of these ideas to the triton ( 3 H) and hclion ( 3 He) was not immediate be¬ 
cause a system of three nucleons collapses under the sole effect of attractive, non-derivative 
contact forces ra. The solution within pionless EFT is the existence of a non-derivative 
three-body interaction mm at LO. This force provides not only saturation, but also a 
three-body parameter which correlates certain observables such as the triton binding energy 
and the doublet neutron-deuteron (nd) scattering length (“Phillips line” [2D]). This frame¬ 
work recovers Efimov’s universal approach to the three-nucleon problem [2X1 [23] . With recent 
progress, pionless EFT also allows for an elegant and efficient fully perturbative treatment 
of contributions beyond LO [TB1 124H2E] . 

An important question is how far up the nuclear chart this EFT applies. Nuclear density 
tends to increase with nucleon number A, implying larger typical nucleon momenta within 
the nucleus. Calculations suggest that pionless EFT holds for A = 4 [EMO] and perhaps 
up to A = 6 systems [301 El] without a four-body force at LO. As the pion mass increases, 
its range of applicability increases, and this framework has been established as a powerful 
tool that can be used to analyze and supplement calculations of light nuclei directly from 
lattice QCD [D211D4] . (See Refs. IPol!D7] for earlier work on pionless EFT for unphysical 
quark masses using input from chiral potentials.) 

The application of EFT to nuclei requires an understanding of the importance of Coulomb 
and other electromagnetic effects. The long-range nature of the Coulomb force implies that 
it becomes dominant at very low energies, i.e., precisely where the EFT is supposed to 
work best. To describe scattering in this regime, one thus has to resum Coulomb effects to 
all orders to recover the Coulomb-modified effective-range expansion. In the proton-proton 
(pp) sector of the pionless theory this was first carried out by Kong and Ravndal [38] [39], 
with subsequent discussions of the charged two-body sector given in Refs. 1401 !4o] . An early 
attempt at pd scattering was made in Ref. [40], and extended to lower center-of-mass energies 
in Ref. [TTf. A calculation of the Coulomb-modified pd scattering length, with the consistent 
use of a Yukawa-screened Coulomb potential in momentum space, has been presented in 
Ref. [25]. 

The 3 H- 3 He binding energy difference has been studied using pionless EFT in a number 
of papers. Systems of three and four nucleons including the Coulomb interaction as part 
of an effective potential were first analyzed by Kirscher et al. using the resonating group 
method [281 l30i 1X8] . Ando and Birse [49] carried out a momentum-space LO calculation 
which was nonperturbative both in the sense that it looked for the bound state as a pole in 
the pd doublet-channel amplitude, as well as in the electromagnetic sector, where Coulomb 
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effects were included through a fully off-shell Coulomb T-matrix, using methods developed 
by Kok et al. in Refs. [50, [51]. A subset of the present authors presented a calculation 
of pd scattering and the bound-state regime in Refs. [13 E3 E3]. Of these, Refs. [13 E2] 
included a perturbative framework using trinucleon wave functions. An updated version of 
this calculation that corrects some issues of the previous approach has recently been given in 
Ref. [50]. That work also includes a nonperturbative treatment, in which all 0(a) Coulomb 
diagrams are resummed to all orders. This calculation was similar to that of Ref. [09] but 
found that the full Coulomb T-matrix is not necessary for an accurate description of the 
3 He bound state. 

At LO in the strong interactions, the perturbative and nonperturbative results of Ref. [50] 
were found to agree with each other (as well as with the experimental value for the 3 H- 3 He 
binding energy difference) to within roughly 30%. While at first sight this seems fine if one 
keeps in mind the expected uncertainty based on the EFT expansion, the relevant parameter 
in this case is actually not the Q/ A^ of the strong sector (with the typical low-momentum 
scale set by the deuteron binding momentum, Q ~ ~ 45 MeV and the breakdown scale 

Af ~ ~ 140 MeV), but rather the aM^/Q scale (with M'n ~ 940 MeV) set by the 

Coulomb contributions. Since in the bound-state regime the momentum scale is set by the 
trinucleon binding momentum, ~ 80 MeV, Coulomb effects are expected to be a small 
perturbation. Based on this, one should expect better agreement between perturbative and 
nonperturbative calculations. 

The purpose of this work is to investigate a rearrangement of pionless EFT that explores 
the role of the increasing nuclear binding momentum in the simplest context, the trinucleon 
systems. We develop a new formulation to include Coulomb and other electromagnetic effects 
in nuclear ground states using perturbation theory. There are two important ingredients to 
this procedure. First, we introduce a new counting scheme that takes as LO the spin-singlet 
channel in the unitarity limit (infinite scattering length) and only includes the finite 1 S'o 
scattering length as a perturbative correction. We will show with the 3 H binding energy 
and with doublet nd scattering phase shifts that deviations from 1 So unitarity are indeed 
small. In the case of pp scattering at very low energies, the scattering-length term is iterated 
along with Coulomb effects so that renormalization can be carried out by matching to the 
Coulomb-modified effective-range expansion. This is essentially what was introduced in 
Ref. [39], but here we isolate the contribution from the divergent single-photon bubble, 
which was missed in the perturbative calculation of Ref. [50]. This new scheme is the second 
ingredient which then allows us to use the counterterm fixed by pp scattering also in the 
perturbative calculation of the 3 H- 3 He binding energy difference at next-to-leading order 
(NLO) in the new counting. The result differs from that of the nonperturbative calculation 
by much less than the 30% previously obtained, and is slightly closer to the experimental 
value of the binding energy difference. At NLO we end up about 1.5% off the 3 He binding 
energy. 

Besides finally obtaining a complete calculation of the energy splitting that includes 
Coulomb effects only as perturbative corrections, we also employ a fully perturbative treat¬ 
ment of corrections in the strong sector. This was already done in Ref. ra. which showed 
that a new, isospin-breaking three-body counterterm is needed for renormalization in the 
presence of range corrections when two-body Coulomb effects are resummed to all orders. A 
further motivation for the current paper is to revisit this issue with a completely perturbative 
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inclusion of those contributions as well. 1 

This paper is structured as follows. In Sec. |TT] we discuss the basic formalism of pionless 
EFT with explicit isospin breaking. The two-body sector and in particular our new method 
to treat the pp channel is described in detail in Sec. Ill In Sec.|IV|we discuss the implications 


for the three-body sector. The new results for the 3 H and 3 He binding energies are presented 
in Sec.[V} We conclude in Sec.[Vl]and provide technical details about the divergent Coulomb- 
bubble diagram in the Appendix. 


II. EFFECTIVE LAGRANGIAN 


We are interested in describing nuclear bound states in terms of the most general dynamics 
consistent with QCD symmetries built from a nucleon held N of mass Mjy. We use a modified 
and extended version of the notation and conventions in Refs. [52, [5HJ . To NLO we write 
the Lagrangian as 

C = jV 1 ^i-Do + ^ ^ 2d ^ 2t > (1) 

where D M + ie A^Q^ includes the direct nucleon coupling to the electromagnetic held 

via the charge operator Qn = (1 + ts)/2. 

It is convenient to express the interaction terms via auxiliary dibaryon helds d l and t A in 
the NN channels projected with 

pj = (X Vr 2 /V8 , P a = ct 2 t 2 t a /V 8, (2) 


respectively, where lowercase (uppercase) letters are spin-1 (isospin-1) indices. While Qn 
ensures that Coulomb photons couple only to protons, the isospin-1 basis used for the 
spin-singlet held t A otherwise mixes contributions from protons and neutrons. If we take 
A = 1, 2, 3 to be an index in the Cartesian basis, then the np configurations are completely 
contained in the A = 3 component, whereas nn and pp are obtained from linear combinations 
l±i2. This corresponds to using a spherical basis A = —1, 0,1, where one would immediately 
have the desired separation. We thus define new projectors 

pp = ( p ? - if, 2 ) . p° = f, 3 . pp = (f, 1 + if, 2 ) (3) 

that correspond to pp , np and nn configurations, respectively. Since the transformation from 
P A to P A is a unitary rotation, the normalization is automatically correct: 

Tr ((PPP?) = \iiB ■ (4) 


In terms of these helds and projectors, the two-body interactions are 


Co d = -d 


ffd, + Cd 




d l + y d [dC (N T P l d N)+ h.c.] 


(5) 


1 The viability of perturbative Coulomb exchange for the trinucleon system was also investigated in inde¬ 
pendent work |.15] . which appeared shortly after submission of our manuscript. 
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D 2 


AM 


N 


t° - t ~ lf 


®t,pp T Cj 


■PP 


i-Dn 


D 2 


-i 


-t +1 1 


,rm “1“ Ct,nn i-Do + 


D 2 


AM, 
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t +1 + y t 


AMjy 

i" (n t P i a n) +h.c. 


( 6 ) 


The covariant derivatives include the appropriate charge operators Q. To keep the notation 
as simple as possible, we use the plain subscript to refer to the np channel from here on 
and use further qualifications only to denote pp and nn (the latter is not explicitly considered 
in the rest of this paper, except at the end of Sec. [Vj). The parameters a d and a t (~) are 
related to the respective scattering lengths; each is actually a sum of contributions from 
various orders: 


a t (,-•) “ 


(o) 

d 

( 0 ) 


r (1) 

’d 

r (l) 


t + + 


(7) 

( 8 ) 


A more detailed discussion is given in Ref. |5T] • Departing from what is used there we follow 
here Ref. [56] and simply set 


2 2 

y d = yt = 


An 

M n 


(9) 


At the same time, the new parameters q and c^..) have been introduced to incorporate 
effective-range corrections, cf. Fig. [2j starting at NLO: 


°d — 4^ ■+— j ( 10 ) 

= c t (1) H-• (11) 


In writing Eqs. ([8]) and ( [If] ) we imposed isospin symmetry at the lowest order of each 
parameter. The reason is that isospin violation comes from either electromagnetism or the 
up-down quark mass splitting. These are associated with mass scales aM N ~ 7 MeV and 
m u — rrid ~ 3 MeV that are small compared to the breakdown scale A^. 

The three-nucleon interaction that is needed already at LO to renormalize the doublet- 
channel amplitude ra can be written as (Ml [57] 


£3 = [yl tftdW + y 2 t t A h B T A r B - y d y t {d l h A aW A + h.c.)] N , (12) 

where the coupling h is also split up in various orders, 


h = /i (0) + h m H-. 


(13) 


As we are going to show below, there is no need for an isospin-breaking three-body force to 
NLO in our expansion. Higher-order terms, including further isospin violation, are briefly 
discussed in Sec. El 


III. TWO-BODY SECTOR 

In this section we use the Lagrangian from the previous section to derive the two-body 
propagators and amplitudes which are basic ingredients of the three-body calculation de¬ 
scribed in the next section. 
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Figure 1: Full dibaryon propagators in (a) the 3 Si state (i.e., the deuteron) and (b) the 1 So state. 
A single solid line represents the propagation of a nucleon, a double dashed line a bare d propagator, 
and a line of circles a bare t propagator. 

A. Spin-triplet propagator 


The dibaryon residual mass ad represents the physics of the triplet NN scattering length 
or, alternatively, the deuteron binding momentum 7For momenta Q ~ 7^ <C Ay, the 
standard power counting of pionless EFT applies [2I5U7U9]. The bare LO dibaryon prop¬ 
agator is simply —i / a^\ and it has to be dressed by nucleon bubbles to all orders in order 
to get the full LO expression. Summing up the geometric series shown diagrammatically in 
Fig- Ufa) gives 

(14) 


iAi 0) (po,p) = 


a 


d 0) + PdMPo, p ) 


with the bubble integral 


r-A 


Iq(PO,P) = M N 


d 3 g 


(27t) 3 M N p 0 - p 2 /4 - q 2 + 1 £ 

Mn I 2A 

7r 


47T 


y - M N p 0 - ie ) + 0(1/K) . (15) 


Here, we have used a simple momentum cutoff A for regularization (not to be confused with 
the physical scale Ay- at which pionless EFT breaks down because new dynamical degrees 
of freedom are resolved). To get expressions in the more commonly used power divergence 
subtraction (PDS) scheme [5], one has to replace 2A/7T —>■ p, where p is a scale introduced 
through dimensional regularization. Of course, in a properly renormalized theory the choice 
of regulator is arbitrary. While power counting is clean and transparent in the PDS scheme, 
we find it convenient here to work with the simple cutoff instead, because this regulator is 
more straightforward when it comes to the consistent treatment of Coulomb contributions. 


For completeness we note that Eq. (15) is valid up to corrections proportional to inverse 


powers of the cutoff, which we neglect here. 

Range corrections are accounted for by the dibaryon kinetic parameter q, which is in¬ 
cluded fully perturbatively here. At NLO, we consider one insertion of the operator shown 
in Fig. 2^a) between A^s. Depending on the renormalization procedure, the insertion of 
c'P might require a concomitant insertion of a^}\ 


iA^GocbP) = iA^(p 0 ,p) x 


(o), 


-\a 


(!) :J1) 


1C 


Po 


P^ 


4Af, 


N 


x CWp) 


(16) 


In the spin-triplet channel, we use the effective range expansion around the deuteron pole 
for renormalization and thus require that 

4tt \ 

- i m = = e/M N , P = 0 ) = ( i7) 


6 
















rsj 


A.Z. 


(-i Cd) 



4 Mjv 


* 


(-ic t ) (p 0 - ife) 


(a) 


(b) 


Figure 2: Dibaryon kinetic-energy corrections in (a) the 3 Si 


with the perturbative expansion 


1 

k cot 5 d (k) — i k 


1 

-7d + y(^ 2 + 7d) +- [k 



1 Pd^+Jl 
2 7 d + ifc 



(18) 


where 7 ^ = \/M N E d ~ 45.7 MeV [58] is the deuteron binding momentum and p d ~ 1.765 
fm [5P] the deuteron effective range. 

At LO, the effective range p d ~ 1/A^ does not contribute and one simply has 



2A 

- Id- 

7 r 


(19) 


The corresponding propagator A^ is given by the term outside the parentheses in Eq. (18), 
up to an additional minus sign. 

The corrections in the parentheses come from subleading orders. The first of these are 
0(Q/ A^y^/A^), or NLO, and 


(1) _ Pdl d (1) _ M N p d 
— 2 ’ C d — 2 


( 20 ) 


B. Spin-singlet propagators 


In the spin-singlet channel, in the absence of Coulomb effects, we can follow the same 
procedure as in the spin-triplet channel (see Figs. 0 b) and | 2 ](b)) , and arrive at the LO 
propagator 


W 0) (po,p) 


—1 

4 °’ +Si 2 -fo(Po.p) 


( 21 ) 


At NLO, we have 


iA t ( 1 ) (p 0 ,p) 


iAi (0) (p 0 ,p) x 





x iA[ 0 ) (p 0 ,p) • 


( 22 ) 


Renormalization is then performed by matching A t to the 1 S'o effective range expansion 
around zero momentum, 


k cot 8Ak) =- 1 — -k 2 + Oik 4 ), (23) 

a t 2 

where a t — —23.7 fm and r t ~ 2.73 fm [60] are the rip scattering length and effective range, 
respectively. As for the triplet, r t ~ 1 /A^, but |a*| l/y^. 
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1. Standard approach 


Typically, one is interested in low momenta Q ~ 1/1 | <C A^ and demands that at 

leading order the scattering length is reproduced. Effective-range and higher corrections are 
included perturbatively, since they are C>((2/A^, l/(|a t |AM). Up to NLO in this standard 
scheme one sets 

(o,st) 2A 1 (i s t) (ist) M N r t 

- -, - 0 > c - —— 

7 t at 2 

There is no adjustment of cr t when the effective-range expansion is performed around the 
zero-energy threshold. 



2. Unitarity limit 


The standard approach makes sense for momenta Q ~ 1 /1 a* |, so at LO we have contri¬ 
butions from both the unitarity cut and the scattering length. However, if we are interested 
in momenta 1 /1 a £ | <C Q <C A^, we are close to the unitarity limit and take instead 


cr 


(o) 

t 


2A 

71 



1 

a t 



M N r t 

2 


(25) 


so the actual finiteness of the scattering length only enters as a perturbative correction. The 
difference with respect to Eq. (24) is a result of \a t \ l/y^. For example, for Q ~ we 
are performing an extra expansion in the ratio l/{\a t \^d) •C 1. Meanwhile, other effective 
range parameters require no special treatment, being similar in both channels—for example, 

r t ~ p d - 

Of course, to the extent that this expansion works one might as well resum ay 1 and use 
Eq. (24) with = aj + a\ ]> over the whole range Q -C A^. However, there are also 

advantages in keeping a strict ordering. It makes clear, for example, that the singlet LO is 
parameter free. Also, as we show shortly, it matches well with the perturbative expansion 
of electromagnetic effects. 


C. Coulomb insertions 

Simple dimensional analysis shows that the Coulomb expansion is in powers of aM^/Q, 
while other electromagnetic corrections are suppressed by at least (Q/ A^) 2 . To NLO we 
need only keep the contribution from static Coulomb photons. In this case, the matching 
discussed above is correct for the np part of the 1 S' 0 dibaryon. If one neglects strong isospin 
breaking, it can also be used to describe the nn component. For pp configurations, however, 
one has to use the Coulomb-modified effective range expansion [U SSj, 

Cl (k cot S ttPP (k) - i k) + aM N H(ri) = —— + ^k 2 -|- (26) 

ac 2 

because electromagnetic effects dominate the very-low-energy scattering regime, Q < aM N . 
These are encoded in the Coulomb parameter 

ctM]\r 

~ 2k 
8 



(27) 







the Gamow factor 


(28) 


and the function 


2 = 27T7/ 

e 2nri _ 


H(v ) = ^O?) + ^ 


log(i?/), 


(29) 


with 0 the derivative of the Euler Gamma function. Here a E — —7.81 fm and rc — 
2.79 fm [61] are the Coulomb-modihed scattering length and effective range, respectively. 


One arrives at Eq. (26) after subtracting the pure Coulomb amplitude with Coulomb phase 
shifts 

exp(2i(r 0 ) = T(1 + h?)/r(l - irj) (30) 

from the full amplitude. 


1. Kong+Ravndal approach 


Since aM N ~ 1 /1 a f |, describing physics at the scale of the 1 S' 0 virtual state requires 
also a resummation of Coulomb exchange. This has been studied in detail by Kong and 
Ravndal in Ref. [39] in a setup without dibaryon fields. We refer to that reference for details 
and note here that the relation of our parameter to the Co of Kong and Ravndal is 

agf = -47T/(M iV Co). 

The key ingredient is the fully dressed Coulomb bubble shown in Fig. [3| From an eval¬ 
uation of the Coulomb Green’s function, Kong and Ravndal find that it is given by the 
divergent bubble integral 


Mk) — Mn 


d 3 g 27 rr](q) 


(27r) 3 e 2nr >^ 


1 k 2 — q 2 + \£ ’ 

and they separate the divergent part with a subtraction at k = 0: 


Jo(k) — M]y 


d 3 g 27r?7(g) 


k 2 


(27t) 3 e 27 "?^) — 1 q 2 [k 2 — q 2 + ie) 


-M, 
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d 3 g 


(27t) 3 — 1 q 2 


= -CM + J 0 

where the finite piece is 


div 


Jo(k) = 


aMJj 
An 


H{r ?) • 


With a simple momentum cutoff, the divergent part is 

2A 


jdiv _ 

J 0 ~ 


Mjy A aMfj 

‘In 2 47T 


log 


aM, 


N 


C E + <D{ 1 /A), 


(31) 


(32) 


(33) 

(34) 


where C E ~ 0.5772 is the Euler-Mascheroni constant, and we neglect higher-order correc¬ 
tions in 1/A. 

Resumming now this dressed bubble along with , the singlet dibaryon propagator 
in the pp channel becomes 


iA 


(0,st) 

t,pp 


(Po,p) 


—1 


AO,st) 

U t,pp 


2A 

n 


T aM N 



2A 

aMjy 



aM N H(rf) 


(35) 
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Figure 3: Fully dressed proton-proton bubble. A wavy line denotes a Coulomb photon exchange. 


with 


1 1 = '/(i v P-/4 - -Vv/'m - i " j 


(36) 


Matching to the Coulomb-modified effective range expansion is performed via the part of 
the T-matrix where Coulomb interferes with the short-range interactions, 


i T sc (k) = C 2 e 2ia °(iy t ) 2 iA ttPP (po = k 2 /M N ,p = 0) = i 


47T 


3 2i (tq 


k cot 8t,pp{k) — i k 


(37) 


The additional factors here compared to the np component, Eq. (17), arise from the inclu¬ 


sion of initial and final-state Coulomb interactions in order to get the amplitude from the 


propagator. Combining this relation with Eq. (26) one finds that both the Garnow factor 


and the pure Coulomb phase shift drop out and one arrives at 


(0,st) 


1 2A 

-1- 

a c 7r 

Range corrections were also considered in Ref. 


a ‘.PP' = - rz + IT - aM * ( l0 S 


2A 


c f 


N 


(38) 


2. Separate leading-order resummation 


Here we develop a new approach that allows us to consider a fully isospin-symmetric 
leading order in the spin-singlet channel, including the pp part. Since ac is still large 
compared to the typical nuclear length scale set by the inverse pion mass (1/m w ~ 1.4 fm), 
in the momentum window aMjy < 1/ac < Q C Ajf we remain close to unitarity in the 
singlet channels. We should be able to treat Coulomb perturbatively along with finite 
scattering-length corrections. The key idea is that with the new counting, the LO singlet 
propagator now behaves exactly like l//c, he., it has the same infrared behavior as Coulomb 
contributions, for which the relevant parameter (in the pp system) is r] = aM^/{2k). 

At LO we thus have, with an isospin-symmetric satisfying Eq. (25), 


iA S P (Po,p) = iAS 0 ) (p 0 ,p) = 


y'r - M N p 0 - ie 


(39) 


At NLO, we need perturbative insertions not only of a\^ p (Fig. |d|a)) and c t V> (Fig. 


but also of a single-photon exchange, see Fig. |5](b). Using the notation of Ref. 
the single-photon piece h/o(fc) and find for the correction to the pp propagator: 


we call 


iA S P (Po, P) 


iA^(p 0 ,p) 
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-icr, 


(i) 

t,pp 


1C 
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p^ 


AMn 

iy 2 5I 0 (iy/p 2 /4 - M N p 0 - \e 


x iAi°pp(po, P) 


(40) 
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(a) 


(b) 


Figure 4: Corrections in the pp channel: (a) Coulomb-corrected scattering length; (b) one-photon 
exchange. 



Figure 5: Dressed proton-proton bubble excluding the empty piece. 


In order to match to the Coulomb-modihed effective range expansion, we work in the pp 
center-of-mass frame in the remainder of this section. From the expression for SIo(k) in 


Appendix Al, we find for p 0 = k 2 /M Nl p = 0, and r] as defined in Eq. ( |27| ) that 
SI 0 {k) = 


aM^ 

An 


log i rj + log 


2A 


cuMn 


~ Cr 


+ 0(1/A ), 


(41) 


where — 1.119 and we again drop terms proportional to inverse powers of A. The log 
divergence is absorbed into a[ p P , whereas the log(i? 7 ) constitutes the one-photon contribution 


to H(rj) as defined in Eq. ( |29| ). 

In order to hnd the renormalization condition for cr^ pp) we now consider smaller Q, which 

requires the resummation of both a^ pp and Coulomb exchange. We determine a^ p from this 
calculation, and then use it in the regime where Coulomb is perturbative, to order a. For 
this, it is important to be consistent about which finite terms get absorbed into along 
with the logarithmic divergence. 

The resummation of Coulomb produces a new dressed bubble, shown in Fig. [5] The im¬ 
portant point is that this new dressed bubble excludes the empty piece (no photon exchange 
inside the bubble) because that has already been resummed at LO. We denote by 8Jo(k) 


the part remaining after subtracting the single-photon piece 5Io(k). In Appendix A4 it is 
demonstrated that 


SJ 0 (k ) = 


aM 2 N 

An 


w*) + ^+ c * 


A, 


47T 


where Ca 


0.579. Hence, we hnd for the new dressed bubble 

2A 


SI 0 (k) + fiJo(k) — 


aMff 

An 


log 


ccMjv 


C c -C A - H( V ) 


+ k. 

An 


(42) 


(43) 


The resummation of a^ pp and the new dressed bubble is now straightforward. We consider 
here explicitly the on-shell case in the center-of-mass frame —po = k 2 /Mpj, p = 0—and hnd 


(A:) = iA<$„(*;) + iA™ (*) 


( 0 ) 


-icr. 


(i) 

t,pp 


i yplo(k) - iyfSMk )) iA® (k) + 



( lo S 


2A 


aMj 


N 


- c c - C A - aM N H{rf) 


(44) 
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where the explicit imaginary part i k cancels. Now, from Eqs. (37) and (26), 


CT- 


(i) 

t,pp 


= -ctMtv ( log 


2A 


ac 




aM, 


- Cr - Ca 


N 


(45) 


Our renormalization of a f ^ p + aj'J p is consistent with Eq. (38), but with a different constant 
finite piece. As discussed in more detail in the Appendix, the reason for this change is that 
we have isolated here the divergent piece—the bubble with a single photon exchange—and 
only regularized that, whereas Ref. [39] regularizes the fully resummed bubble—including 
the empty part—at once. Effectively, this amounts to using different regularization schemes. 
Our approach has the advantage that it allows for a consistent matching between the per¬ 
turbative and nonperturbative regimes. 

The expression for the pp phase shift in the regime aM N < l/|ac| <C Q <C can now 


be found by inserting Eqs. (39) and (40) into Eq. (37): 


k cot S t ,pp(k) = —— + aM N CA + 77 k 2 + aM N log ( 
ac 2 \ 


aM, 


N 


2k 


+ 


(46) 


This is consistent with a direct expansion of Eq. (26) in powers of aM^/Q and l/{\ac\Q)- 
Using 


= ~ ~ C 'E + 0(rj ), 
V 


one finds Eq. (46) with Ca = Ce and 


r NLO _ 

r c — r t ■ 


(47) 


(48) 


Moreover, the same result as in Eq. (46) can be obtained in a direct calculation that includes 


all O(a) diagrams contributing to pp scattering, treated in perturbation theory. 

As discussed in Appendix A4 a numerical calculation indeed yields Ca ~ 0.579, very 
close to Ce ~ 0.5772. Equation (46) gives circumstantial evidence that in fact Ca = Ce, 


and we expect that increased numerical accuracy would reveal the same. Lacking a formal 
proof, however, we keep Ca in the expressions here and merely note that the deviation from 
Ce is negligible for the results presented later. 


Equation (46) is in agreement with the Nijmegen phase-shift analysis up to about 
15% for k > 60 MeV. In particular, it captures the correct slope at intermediate momenta, 


which is a consequence of the fact that Eq. (48) works at the 3% level. The assumption of 
isospin symmetry in c^, which means that we can neglect the splitting in the spin-singlet 
ranges, is a good one. 


IV. THREE-BODY SECTOR 

We now consider the three-body sector in order to calculate the 3 H and 3 He binding 
energies. These bound states arise in the nd and pd spin-doublet channels, respectively. Since 
we do not reorganize the pionless EFT expansion in the NN spin-triplet channel, no changes 
are needed in the nd and pd spin-quartet channels with respect to Refs. 

We take the point of view that the scale that characterizes these bound states, the triton 
binding momentum 7 t, is comparable to the deuteron binding momentum 7 d, but both are 
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much larger than qlMn , l/|cit|, and l/\ac\- Thus the bound-state energies can be expanded 
not only in powers of Q/A^ but also of K 0 /Q, where K 0 ~ aM N ~ l/|a f | ~ l/|ac|- 2 For 
simplicity we pair the two expansions by taking Q r\j (*w ) 1/2 - 


More formalism 


According to the standard power counting PEUIH [12] one-nucleon exchange between 
nucleon and dibaryon has to be treated exactly. To NLO it is sufficient to consider the 
S'-wave projected one-nucleon-exchange diagram at energy E, 


TS (V , A _ 1 ^ fk 2 + P 2 ~ M n E - ie 
K s [E; k,p) = — Q o - - - 

kp V kp 


2 kp 


log 


k 2 + p 2 + kp — M n E — ie 
k 2 + p 2 — kp — M n E — ie 


(49) 


where k (p) is the incoming (outgoing) center-of-mass momentum. It is well known 
that at LO the resummation of one-nucleon exchange is renormalized by the three-body 


force given in Eq. (12), with 


h^ = 


M n H{ A) 
A2 


(50) 


where A is now the momentum cutoff applied in the three-body equations discussed be¬ 
low, and H( A) a known log-periodic function of the cutoff that depends on a three-body 
parameter A*. Here we follow the procedure employed in Ref. [3] and much of the subse¬ 
quent literature, in which the two-body cutoff is taken to be very large and 0(1/A) terms 
in Eqs. (15) and (41) are neglected. 3 


In order to calculate bound-state energies and perturbative corrections to them, we use 
the formalism and notation of Refs. [52, [54] to study the Nd doublet channel. We introduce 
a three-component vector of vertex functions in channel space, 


B s = (H s d ’ a , £> d,bl , H s d ’ b2 ) 


(51) 


where H d,a corresponds to the deuteron channel and H d,bl and £> d,b2 are the np and pp/nn 
components of the 1 S'o multiplet, respectively. These are given by (properly normalized [47]) 
solutions of the homogeneous equation 


B s = (. KD) <g> B s 


E = —E 


B 


(52) 


where <g) represents an integral over the intermediate momentum, D is a diagonal matrix of 
dibaryon propagators written in the form 


D d ,t( E i q) = A d,t [e - q'j , (53) 


2 Note that Ho is an extension of the original definition HIS] to include additional scales, in particular 
aMjv. 

3 This is effectively equivalent to using dimensional regularization in order to renormalize the two-body 
sector first. 
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and 


K = 


«.(*■•+ ^») 
\9«(2K.+ t §P) 


-9u(K s + ^) 

-9u( 2 K. + i%p) - a „i§^ , 


(54) 


The factors gdd etc. contain the coupling constants yd and y t . In our present conventions, 
we simply have 

gdd = gtt = 9dt = 27T. (55) 


For illustration, we note that with a single channel and no three-body force, the homogeneous 
integral equation written out explicitly would be 


B s (p) 


1 
7r 


dqq 2 K s (E',p,q) A E 


q 2 . \ 

2M n ’ 7 


5 s (g), 


(56) 


where A denotes a generic dibaryon propagator, and A is the momentum cutoff employed 
in the three-body sector. 

To study doublet-channel nd scattering, we have to solve the Lippmann-Schwinger equa¬ 
tion. This can be done in a simpler two-channel formalism, where it becomes 


9dd (K. + *%£) \ 

-9dt (3/Ts + ) 

(-g dd Dd (a; + g dt D t (dK s + \ 

+ (3 K s + (/F s + J ^ U d ’ b J ' 

Again, the notation is the same as in Refs. [52, !54j. The nd phase shift is then obtained 
from the on-shell amplitude, 



5n-d(k) = i log[l + 2 -^^Z,T s d '\E k - k, k) 


Ei. — 


3 k 2 


Id 


4Mjv Mtv 

where Z Q is the deuteron wavefunction renormalization determined by 

d 1 


^o' 1 = i 


dp 0 A d {p) 


Po=-ivA'P=0 


The scattering length is determined by the on-shell amplitude in the limit k —> 0, 

o Mn 


@"n—d 


37T 


lim^ZoT^Ef.k.k) 


(58) 


(59) 


(60) 


B. Perturbative range corrections 


Before we discuss Coulomb matrix elements between vertex functions, it is instructive to 
consider effective-range corrections in this framework. As we emphasized in Sec. 


Ill 


range 
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Figure 6: Range corrections contributing to the trinucleon binding energy in perturbation theory. 
A shaded oval represents a vertex function. 


corrections are 0(Q/ A^). In analogy to Eq. (18), we write 
D d {E-q) = Df{E- q ) + D { p{E- q ) + ••• 


-1 


-Id + v 73 ? 2 / 4 - M nE - 


x 


1£ 


x + Pd (3g 2 /4 - M n E - yj) + 


-Id + \/3g 2 /4 - M n E - ie 


, ( 61 ) 


with an analogous expression for the spin-singlet part. Suppose now we have solved the 
homogeneous equation ( [52] ) at leading order. The binding-energy shift due to the deuteron 
effective range, shown in Fig. [6]a), is then given by 


= d 11 2 l B “’"(9)| 


r A 


(3g 2 /4-M^- 7 g) 

_\ 2 

—'id + y/ ; WV4 - MnE - 1 : j 


(62) 


More details about this can be found in Refs. [2TJ |26]. In particular, we note that the 
correction h ^ to the three-body force is fitted to keep whatever physical parameter has 
been used at LO (triton binding energy or doublet-channel scattering length) unchanged. 

In the standard counting used previously (finite scattering lengths at LO), there is a 
contribution analogous to Eq. (62) with y^ —>■ 1 /a t in the denominator (and no in the 
numerator). In our new scheme with the spin-singlet LO in the unitarity limit, the range 
corrections in the spin-singlet channels are particularly simple. For example, from Fig. [6](b) 
we have 

AD 3r t M 


A E, 


r t ,bl 


= 4P / d« 2 |A bl (9)|' 


(63) 


and an analogous expression that involves NLO corrections that are linear in the range 

have been known for a long time to generate cutoff dependence that can be compensated by 
h D pJJ mm- F° r three bosons at unitarity, this divergence was discussed in Refs. [631IM] , 
At the same time, we have of course now perturbative insertions of the scattering length, 
e.g., 

3 r '. , | e ?- bl ( 9)| 2 




3? 2 /4 + M n E 


(64) 


B 


which are corrections of 0{^q/Q) relative to LO. 
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(a) 



(b) 



Figure 7: Convergent diagrams contributing to the 3 H- 3 He binding energy difference in perturba¬ 
tion theory. 




Figure 8: Divergent Coulomb bubble diagram (a) contributing to the 3 H- 3 He binding energy 
difference along with the associated counterterm (b). 


C. Coulomb matrix elements 


Range corrections as discussed in the previous section apply in general to both nd and 
pd systems. In the bound-state regime, they simply correspond to matrix elements between 
trinucleon wavefunctions that are diagonal in momentum space (only one loop integral is 
required to calculate them) as well as in cluster-configuration (channel) space [56]. 

Now we want to include Coulomb corrections, which are 0(Hq/Q). In general, contribu¬ 
tions to the 3 H- 3 He energy splitting A E can be non-diagonal in both spaces. We get such 
contributions when we calculate A E in perturbation theory, which comes from diagram 
topologies shown in Figs. [7] and [8} This approach starts by taking the trinucleon state to be 
the triton in a three-channel formalism. This way, one can easily isolate the channel that 
corresponds to the pp configuration in 3 He. Such a calculation was carried out in Ref. (|53J for 
the diagrams in Fig. [7J which give convergent results. We briefly summarize this calculation 
here before getting to the new diagrams in Fig. [8j 

The diagram shown in Fig. [T](a) is still diagonal in cluster-configuration space. In order 
to calculate it, we need the kernel function corresponding to the photon being exchanged 
between the deuteron bubble and the individual proton. It is given by m 


-ftbubbi e{E]k,p) =—aM N [ dcosd t ( ubblc ^. 9 ’ —, k • p = kp cos 9 , (65) 


'-i 


(k - p ) 2 + A 2 


with 


arctan 


-^•bubble(-^O k , p) 


2p—k 2 —k p 

v /3k 2 -4M JV E-ie v /(k-p) 2 


+ arctan 


2k—p 2 —k p 

y/ 3p 2 —4Mj\r E—ie yj (k-p) 2 


V( k “P) 


( 66 ) 
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and where A is a photon mass introduced for regularization in the infrared. In practice, we do 
not perform the angular integral numerically but rather use the explicit S'-wave projection 
given in Ref. [21)] . The contribution to the energy shift is then given by 

A£j§I) = T £ d ?1 q\ £ d q 2 ql Bf‘(q 2 ) Df(-E B , „) 

x AT bubbl ,(£; q u q 2 ) Df(~E B , q 2 ) B}'(q 2 ). (67) 

The analogous diagram with np singlet propagators (not shown explicitly in Fig. [7]) is given 
by essentially the same expression with the replacements Bf ,a —>• Bf ,hl and Dd —>■ D t . For the 
“box” and “triangle” contributions, Figs. [7j(b) and (c), the corresponding kernel functions 
are [54i 


K hox (E] k,p ) = -aM n 

( arctan 


1 

x - 
2 


and 


dcos # 


'-i 


V 


2p 2 —k 2 —kp 


v /3k 2 -4M JV £-ie v /(k-p) 2 


+ arctan 


2k 2 —p 2 —kr 


A /3p 2 -4M Ar E-ie^/ (k-p) 2 


K< 


(out) 


tri 


{E\ k,p ) = —aMjy x - 


K™(E-,k,p) = Ki;r‘(E-,p,k), 


(k 2 + p 2 + k • p - M n E - ie) ^/(k - p) ; 

_A_ 

(k 2 + p 2 + k ■ p — M a iE — ie) 2 

l t n(E; k,p) 


0( A 2 




/ dcos#- 0 ,, „ 

!_ 1 k- + p 2 + k • p — M n E — ie 

(out) , 


( 68 ) 


(69a) 

(69b) 


where the superscripts “out” and “in” indicate whether the Coulomb-photon exchange is 
on the left or right side of the diagram. This notation is taken over from Ref. [53], which 


allowed for incoming and outgoing pd states. The loop function appearing in Eq. (69a) is 
given by 


ltri(E; k, p) = 
X ) log 


2A/k 2 /4 + k p + p 2 
i(k 2 /2 - k • p - p 2 - A 2 - M n E - ie) 


v/k 2 /4 + k • p + p 2 

i(k 2 + p 2 + k • p — A 2 — M n E — ie) 


- log 


+ 2 a/ A 2 + 3k 2 /4 - M n E - ie 
+ 2A 


\J k 2 /4 + k • p + p 2 


(70) 


As for K huhh \ e (E] k,p), explicit S'-wave projections where the integral over cos# has been 
carried out analytically can be found in Ref. [2b|. Resulting contributions to the energy shift 
are of the form 


x K boK (E’, qi ,q 2 )Df\-E B ,q 2 )B?-‘(q 2 ) (71) 
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and 


A - £ S’) = jf d,1? ‘ / A d « 2 fe 2B . d ' M ta) Dl°\-E B , qi ) 

X K[T\E-,<H,q 2 )Df(-E B ,q 2 )Bt'‘(q 2 ), (72) 

with analogous expressions for equivalent topologies but different combinations of dibaryon 
propagators and vertex functions (see Ref. m for details). 

The above summarizes the diagrams included in the perturbative calculation of Ref. 0; 
all these contributions are convergent as the cutoff A is increased. As mentioned in the 
introduction, the contribution from the diagram shown in Fig. |8](a) has not been included 
so far. This diagram is logarithmically divergent, but this is precisely the same divergence 
of the one-photon bubble that we isolated in the new treatment of the two-body sector 

Hence, it can be renormalized by including it together with the 


discussed in Sec. Ill C 2 


counterterm diagram shown in Fig. 8j(b) , which is proportional to a ^ as given in Eq. (45). 
The resulting contribution to the energy shift, written out explicitly, is 


a pd) 

A i|a+b) 


r A 


= ^ L dqq 


K h2 (q)\ 


3g 2 /4 + M n E 


x 


B 


- oiMjy 

ac 


C A + log 


aMjsr 


2 a/ M n E b + 3 q 2 /4 


(73) 


Note that the constant C( drops out here against the same contribution from the photon 
bubble, cf. Eq. (41). Expanding the renormalized pp propagator of Refs. jUl IS] 


m a 


gives Eq. (73) with Ca —> Ce , as expected from their similar values (cf. the discussion in 
Sec. |III C 2h. 


We stress here that the new approach takes all spin-singlet propagators in the unitarity 
limit at LO. In the pp channel, the finite scattering length ac is included together with the 


single-photon bubble contribution, resulting in Eq. (73). At the same time, we also include 


linear insertions of l/a t in the np spin-singlet channel, as given in Eq. (64). 


V. RESULTS AND DISCUSSION 

We summarize our new expansion as follows: at leading order, we include 

• the standard NN spin-triplet (pionless) amplitude (parameter 7 d ), 

• the unitary NN spin-singlet amplitude (parameter-free), 

• a contact three-body force (parameter A*). 

Our new NLO includes 4 

• the effective range in the NN spin-triplet channel (parameter p ^), 

4 These NLO contributions induce corrections to the spin-triplet two- and three-body force parameters that 
already appeared at LO, but these corrections introduce no new parameters. 
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Parameter Value Ref. 


Id 

45.7 MeV 

m 

Pd 

1.765 fm 


at 

-23.714 fm [60] 

n 

2.73 fm 

m 

ac 

-7.8063 fm [61] 


Table I: Parameters used for the numerical calculation. 

• the isospin-symmetric range in the NN spin-singlet channel (parameter r t ), 

• a scattering-length correction to unitarity in the NN spin-singlet np and mi channels 
(parameter a t ), 

• a scattering-length correction to unitarity in the NN pp channel (parameter ac), 

• one-photon exchange (parameter a = 1/137). 

The two-body parameters we use in our numerical calculation are summarized in Table |T} 
For the nucleon mass we take Mn = 938.918 MeV. 

Unitarity in the NN spin singlet at LO means that our results for the binding energies 
and scattering of nuclei differ from previous calculations, for example the 3 H binding energy 
and nd scattering in the doublet channel HM8]. In order to facilitate the comparison with 
existing, standard-LO results, we also consider below an “incomplete” new NLO, where we 
set p d = r t = 0. 

In the spin singlet we perform an extra K 0 /Q expansion on top of the standard Qj 
expansion. We first show that the K 0 /Q expansion gives results that are in good agreement 
with the standard leading-order NN spin-singlet amplitude in the absence of Coulomb effects. 
As A* is varied with fixed NN input, doublet-channel observables change in a correlated 
way. The simplest example is the Phillips line [20.J in the plane of 3 H binding energy and nd 
scattering length, see Fig. [9] Five curves are shown for a three-body cutoff A = 2.4 GeV; 
effects from further increasing the cutoff are negligible. In three of the curves the ranges 
are set to zero. We see that the new LO curve (with a t —> oo) is within 1% of the standard 
LO curve (with a t at its physical value), and the new LO+(incomplete)NLO (with a t —* oo 
at LO and 1 /a t at its physical value treated in first-order perturbation theory, but zero 
range) is closer still. The inset magnifies a region of the plot to show the small differences 
among these curves. This agreement is not fortuitous and survives the inclusion of range 
corrections, displayed in the other two curves. Again, the new LO+NLO curve ( a t —> oo 
at LO, and both physical 1 /at and ranges in first-order perturbation theory) is very close 
to the standard LO+NLO (a* at its physical value at LO, ranges in first-order perturbation 
theory). 

In Fig. [9] we also indicate the experimental values of the 3 H binding energy and doublet 
nd scattering length by, respectively, horizontal and vertical lines. Leading-order curves 
(new as well as standard) lie close to the experimental point. Next-order curves (new as 
well as standard) are small shifts in the direction of data, overshooting a bit. 

We can use either the binding energy or the doublet-channel scattering length to deter¬ 
mine A*, and then the other is a prediction that nearly agrees with data. Here we use the 
3 H binding energy, 77 b( 3 H) = 8.48 MeV, as input. At LO this is done by adjusting M°); at 
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Figure 9: Correlation (Phillips line) between 3 H binding energy (in MeV) and doublet nd scattering 
length (in fm) at LO and NLO. The light (green) and dark (blue) solid lines are the results of the 
standard expansion at LO and LO+NLO. The light (green) dotted line is the new LO, the light 
(green) dashed line is the new LO+NLO with ranges set to zero, and the dark (blue) dashed line 
is the full new LO+NLO. All curves are for a three-body cutoff A = 2.4 GeV. Horizontal and 
vertical (black) dotted lines indicate experimental values for binding energy and scattering length, 
respectively. 


NLO, hS l) ensures that the 3 H binding energy remains at its experimental value. The same 
procedure is used in the standard expansion with a t in the LO propagators, and our values 
for h^ 0 ' 1 " 1 come out very close to results in that approach [26] . Then the nd scattering length 


converges as the cutoff A increases, as shown in Fig. 10, where five curves analogous to those 
in Fig.[9]are displayed. More generally, Fig. [TT] shows the predictions for the doublet-channel 
nd phase shifts at low momenta. Again, the effects of treating the finite value of the singlet 
scattering length in perturbation theory are small. 

Thus, the /Q expansion works quite well in the absence of Coulomb interactions. 
In fact, range corrections seem larger than those from the finite singlet scattering length, 
which suggests that the K 0 /Q expansion works better than the Q/Af expansion. With the 
3 H channel properly renormalized and the three-body force fixed, we now consider Coulomb 
corrections to the Tie binding energy. Because our LO Lagrangian is isospin-symmetric, 
the 3 H- 3 He binding energy difference vanishes in our new LO, but it is a prediction—a 
low-energy theorem—at NLO. 

To gauge the effects of perturbative Coulomb corrections, we show in Fig. 12 the result 
of the new calculation presented here at NLO, as a function of the cutoff. The photon 
mass A has been extrapolated to zero in the same way as in Ref. |5T], that is, a linear 
extrapolation based on the range A = 0.4... 0.6 MeV. All Coulomb effects, including those 
in the pp sector, are included fully perturbatively here, meaning that we only consider matrix 
elements between trinucleon wavefunctions that involve a single Coulomb-photon exchange. 
We find that the inclusion of the renormalized Coulomb-bubble diagram, Fig. [8j ensures 
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Figure 10: nd doublet-channel scattering length (in fm) as a function of the cutoff A (in MeV). 
Numerically, the limit in Eq. (60) has been taken by setting k = 0.01 MeV. Notation as in Fig. [9| 



Figure 11: nd spin-doublet phase shift (in degrees) at LO and NLO as function of the center-of-mass 
momentum (in MeV). Notation as in Fig. |9j 


proper renormalization of the three-body energy. This is in contrast with Ref. [26.], which 
resums some Coulomb contributions already at LO and finds a logarithmic divergence at 
NLO when rc = r t . 

The new contribution also provides a sizable (as compared to the total energy splitting) 
modification of the incomplete perturbative results of Ref. [5l] . It brings the full perturba¬ 
tive result very close to the non-perturbative leading-order calculation of Ref. [54], which 
extended the results of Ref. [TS] to much larger cutoff values. This establishes that Coulomb 
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Figure 12: Results for the 3 He binding energy (in MeV) as a function of the cutoff (in MeV). The 
upper solid curve shows the result of the nonperturbative calculation presented in Ref. [53]; the 
result of Ref. (35] is shown as a red triangle. The incomplete perturbative result of Ref. [S3] is 
given by the lower solid line. Crosses represent the new complete perturbative calculation up to 
linear order in the spin-singlet scattering lengths. For comparison, dotted horizontal lines indicate 
the experimental values for the 3 H and 3 He binding energies. The photon mass A has been linearly 
extrapolated to zero. 


effects really are a completely perturbative correction in the 3 He bound state compared to 
the 3 H. We stress that “leading-order” means something different in that paper than in the 
new approach presented here: Ref. [SI], as pionless calculations preceding it, resums certain 
Coulomb effects and includes the scattering length in the leading-order spin-singlet propa¬ 
gators, whereas we take those in the unitarity limit and only include the finite scattering 
lengths as perturbative corrections. At the same time, our calculation thus shows that this 
is a good approximation, as expected from the K 0 /Q expansion, which embodies the fact 
that these scattering lengths are very large. 

At N 2 LO and higher, isospin-breaking effects from the quark masses and higher-order 
electromagnetic effects will contribute. For example, the 0(a 2 ) Coulomb contribution is 
attractive and tends to reduce the splitting found at NLO. Unfortunately, further shorter- 
range interactions will appear. NN input can, to some extent, be determined from NN 
data. Even in this case, however, there may be an isospin-breaking three-body force needed 
for proper renormalization. At that point, one can no longer predict the binding energy 
difference, unless one can determine this force’s parameter from another isospin-violating 
observable. 

One of the higher-order effects comes from isospin violation in the effective ranges. At our 
NLO there is no isospin breaking in these, so r* = rc and there are no effective-range effects 
in the binding energy difference. This is a feature by construction in our new approach: 
our LO state is isospin-symmetric in the spin-singlet channels, so perturbative corrections 
from isospin-symmetric ranges exactly cancel via the NLO adjustment of the existing three- 
nucleon force that keeps the triton in the right place. Once one includes isospin breaking 
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(b) 


Figure 13: Combined “Coulomb + range” corrections contributing to the trinucleon binding energy 
in perturbation theory, a higher-order effect in our expansion. 


in the spin-singlet ranges, r* 7 ^ rc , at some higher order, one recovers again the linear 
divergence (as a function of the ultraviolet cutoff A) that has been identified in Ref. [26J. 

Further contributions that are proportional to the effective ranges come from the direct 
coupling of photons to the dibaryon fields, which are generated by the covariant derivatives 
in Eqs. ([5]) and ([ 6 ]). The corresponding diagrams are shown in Fig. 13 The expressions 


for these diverge logarithmically as a function of the momentum cutoff A, as identified in 
Refs. [2B, [521 El]- Essentially, the scaling of the diagrams in Fig. 13 is the same as for the 
proton bubble with a single photon exchange that we discuss in Sec. Ill C Whereas in that 
case we had momentum-independent vertices and a single nucleon propagator ~ q~ 2 left 
in each loop (after carrying out the energy integrals), we now get a factor 1/q each from 
the ultraviolet behavior of the dibaryon propagators and trinucleon vertices, respectively. 
Compared to the Coulomb correction (67) the diagrams shown in Fig. 13 are suppressed by 


Q/Ayt- This relative ordering is in fact exactly the same as in previous calculations 
33], but in our new counting scheme it means that these diagrams are N 2 LO. 

Thus, these two divergences associated with the effective ranges appear at higher orders 
in our approach. Note that Ref. [26] contains a third source of divergence which is linear in 
the effective ranges: the interference between non-perturbative Coulomb and perturbative 
range effects. This additional logarithmic divergence occurs because Ref. [26] employs the full 
Coulomb-dressed dibaryon propagator at leading order. We emphasize that this divergence is 
absent in our new approach where all Coulomb contributions are treated perturbatively. We 
see here an example of the more general fact that, when singular interactions are involved, 
the cutoff dependence of perturbative diagrams is not necessarily the same as that of the 
resummed series. As a consequence, as noted above, no new three-body interaction is needed 
for renormalization at NLO in our new counting. 

The discussion here can be generalized to other isospin-violating effects. Our expansion 
for explicit electromagnetic effects is in powers of aM^/Q, which we are counting as Ko IQ 
and pairing, for simplicity, with the standard pionless EFT expansion Q/A 7. In addition 
to photon exchange, there are “indirect” electromagnetic effects that take place at short 
distances and appear in the Lagrangian as interactions among nucleons. Moreover, the up- 
down quark mass difference also generates isospin-breaking interactions. The form of the 
interactions are dependent on the way isospin is broken [551 [6BJ. Since the quark masses 
break charge symmetry (a rotation of 7 r around the second axis in isospin space), their low- 
energy footprints will break charge symmetry as well, at least in first-order perturbation 
theory. In contrast, electromagnetic interactions break isospin more generally. 

The most obvious consequence of isospin-breaking interactions is the neutron-proton mass 
splitting 5Mn. We can estimate this as 5Mn = 0(o.Mn /( 47 r), — m^), where we included 

a 47 t expected from a photon loop. It is well known that these effects have opposite signs and 
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are comparable in magnitude, but, as indicated by this estimate, the quark-mass contribution 
is somewhat larger and makes the neutron heavier. One might worry that this splitting, 
appearing as a mass term in the Lagrangian, should be compared to the kinetic terms shown 
in Eq. ([Tj). However, the nucleon mass splitting term can be removed by a redefinition of 
the nucleon field and once this is done the splitting appears only in the kinetic term 
itself. It is thus 0(5Mn/Mm) relative to leading order, a very small effect. 

The dominant isospin-breaking effects are expected to appear in the short-range two- 
nucleon interactions, represented in Eq. ([TJ) via dibaryon fields. With the choice in Eq. (|9]), 
one counts the dibaryon residual masses as small scales, e.g., a t = C>(Ko)- T° first order 
in the isospin-breaking parameters, a t) nn ~ &t is proportional to m u — m* while a t)PP — a t 
is proportional to both (m u — m d ) and aM N . The largest sizes they are expected to have 
are ( a t>nn ~ cr t )/cr t = 0{{m u - m d )/N 0 ) and (a t>pp - a t )/a t = O(aM N /^ 0 , {m u - m d )/N 0 ). 
The electromagnetic contribution in the pp channel of C?(aMjv/No) = 0(1) is just the one 
required to renormalize Coulomb treated as NLO; the aMjy appears explicitly in Eq. (45). 
This counting is consistent since it yields at/ac — 1 = 0(1), while empirical values give 
a t /a c - 1 ~ 2. 

How we count the quark-mass effects is a matter of choice, since m u — md is at the 
QCD level an independent parameter. The estimate above suggests a t /at, nn — 1 = (D((m u — 
m d )/No) ~ 0.3, again consistent with the standard value a t , nn — —18.7 fm [68], which gives 
— 1 ~ 0.25. However, there are significant uncertainties in the value of a t}nn . For 
example, a value a t)nn — —16.1 fm has also been obtained [62], which would mean more 
significant quark mass effects. Conversely, a value closer to a t would more clearly indicate 
m u — m d as a separate scale, much smaller than No- In Ref. [413] a pionless EFT analysis of the 
trinucleon energy splitting A E was carried out at LO in the standard power counting with 
an additional quark-mass, isospin-breaking NN interaction. In this case, A E is correlated 
with at,nni and one can use the experimental value of the former to determine the latter. 
It was found that a t} nn — —(22.9 ± 4.1) fm, which gives a t /a t , n n — 1 ranging from —0.1 to 
0.25. This suggests that the pairing m u — m d ~ aM N f47 T that one could infer from the 
nucleon mass splitting can be applied to NN interactions as well, which results in quark- 
mass effects about an order of magnitude below electromagnetic ones. We therefore took in 
this paper the standpoint that these are N 2 LO effects and do not contribute to the order we 
were working. Once the most important quark mass contribution is relegated to an order 
where new, undetermined counterterms appear, it is no longer possible to constrain at )nn 
from A E [53 ]. 

For higher terms in the expansion, similar arguments can be used, but now tak¬ 

ing into account that their parameters are determined by the high scales as given by the 
pionless EFT power counting. For example a u m i, c t = (D(M n / A^), so we expect 
( Ct,nn - Ct)/ct = 0((m u - m d )/ A^) and (c ttPP - c t )/c t = 0(aM N / A^, (■ m u - ?n d )/A jt ). These 
relations indicate a suppression of C?(N 0 /A^) relative to breaking in a t (,..). They imply 
rc/rt — 1 = 0(aMiy/Ajt) ~ 0.05, in agreement with rc/rt — 1 ~ 0.02 from empirical values. 
Thus, it is consistent to take the electromagnetic isospin-breaking range as an N 2 LO effect, 
with the prediction (48) valid up to about 5%. As we pointed out above, a linear divergence 
appears in the three-nucleon system which then requires a new, isospin-breaking three-body 
force at the same order. 
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Our result for the 3 He binding energy is 


£ B ( 3 He) LO+NLO 


= £ s ( 3 H) + A E nlo 

= 8.48 MeV - (0.86 ± 0.17) MeV = (7.62 ± 0.17) MeV, (74) 


where we estimated the error in the energy difference as 0(aM]y/Q, (md — m u )/Mo) ~ 20%— 
slightly larger than the ratio NLO/LO. Eq. (74) represents 98.7% of the observed value, 
EB( 3 He) exp ~ 7.72 MeV. There is room for higher-order contributions, but most of the 
splitting is accounted for at NLO. 


We emphasize that the error in Eq. (74) should be understood as a rough estimate of 


higher-order contributions. The above value comes from taking the mean of ac and a t for 
1/Mo. By using the physical value for ac as NLO input, we are overestimating the magnitude 
of electromagnetic effects, because N 2 LO corrections of both electromagnetic and quark- 
mass origin contribute to the pp scattering length. Potential-model calculations of photon 
exchange (for example, Ref. [70j) tend to give a A E consistent with an almost model- 
independent estimate of about 680 keV im but not the full isospin violation in the NN 
scattering lengths. Our result cannot be directly compared with these older perturbative- 
photon calculations because we include also shorter-range interactions—as required to ensure 
renormalization—that give the measured pp scattering length. However, the EFT allows us 
to directly study effects when ac is closer to a t , effectively simulating a potential-model 
calculation with an interaction that does not give the physical splitting in the scattering 
lengths. For example we find that ac = —9.0 (—10.0) fm leads to a trinucleon splitting of 
about 700 (600) keV, and we can thus confirm that the older calculations are consistent with 
an uncertainty of about 20% in ac, in line with our estimate of higher-order contributions 
above. 


VI. SUMMARY AND OUTLOOK 

In this work, we have established a rearrangement of the perturbative expansion in pion- 
less effective field theory that takes the spin-singlet nucleon-nucleon channels in the unitarity 
limit. Not only does this allow us to demonstrate quantitatively that nature is very close 
to this scenario, it also enables us to include Coulomb corrections perturbatively on the 
same footing. An important ingredient to this is the consistent isolation of the divergent 
one-photon piece in the pp sector, which then guarantees the renormalization of the corre¬ 
sponding contribution in the three-nucleon sector. 

By studying the trinucleon bound-state sector, we confirm that the new perturbative 
expansion works very well: both the Phillips line and doublet-channel nd phase shifts are 
barely changed by perturbative corrections that include the actual finiteness of these quanti¬ 
ties, proving that the new leading order is a good starting point for a perturbative expansion. 
By comparison with previous nonperturbative calculations of the energy splitting in the 3 H- 
3 He iso-doublet, we furthermore show that the Coulomb interaction is indeed a completely 
perturbative effect in these bound states. Our new NLO with effective-range corrections 
set to zero is almost on top of previous leading-order calculations that resum certain (but 
not all) Coulomb contributions. The same holds when isospin-symmetric ranges are in¬ 
cluded. There is no new divergence at this order, which allows us to predict the energy 
splitting as AE NLO = —(0.86 ± 0.17) MeV, to be compared with the experimental value 
A£ exp ~ -0.764 MeV. 
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The main point of our reorganization is that, by treating small quantities in perturbation 
theory, we remove inessential parameters from the lowest orders. For example, a 30% accu¬ 
racy for nuclear observables susceptible to pionless EFT should be obtained from only two 
parameters at LO, with four other parameters at NLO ensuring a 10% accuracy or better. 
An investigation of this reorganization in heavier systems is our next goal. 

Further improvement comes from higher orders. Isospin-symmetric effective ranges do 
not contribute to the 3 H- 3 He energy difference. The inclusion of isospin breaking in these 
effective ranges will recover the previously identified divergence generated by this effect; the 
same is true for contributions that are of the same order in a as other Coulomb effects 
included here, but suppressed further by Pd,r t ~ 1/A^. Confirming that these contributions 
are indeed best accounted for at N 2 LO, and carrying out a fully consistent renormalization 
at this order, is beyond the scope of this work and relegated to future investigations. 
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Appendix A: The Coulomb bubble integrals 


1. Single-photon bubble 


Let us consider the bubble integral with a single Coulomb-photon insertion. The original 
integral is given by 


SI 0 (k) = 4t iM 2 n 


' d 3 gi f d 3 q 2 


1 


(27t) 3 J (27t) 3 k 2 — qf + ie q 2 + A 2 k 2 — (qi + q 2 ) 2 + ie 


(Al) 


This has been calculated by Kong and Ravndal 


ever, their result is not entirely correct, so we repeat the calculation in Sec. Ala 


using dimensional regularization. How- 

In 


Sec. Alb we carry out the calculation with a simple momentum cutoff, as used in the main 
part of the paper. 
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a. Dimensional regularization 


Following Kong and Ravndal, we set k 2 = —y 2 . Going to d = 3 — e dimensions, and 
introducing the renormalization scale // and a Feynman parameter, we get 


6 h = ^aMlD) 2 ‘f^^— I'dxl^L -±- 

' 2 / J {2n) d ql + A 2 J Q J (2n) d [ x (^ 2 + qj) + (1 - x)("f 2 + {(h + <h) 2 ) 


d d qi 


= AiraMlj (^ 


2 e 


r h + D r' dx t id i 


( 4 vr)^r( 2 ) Jo ' J ( 2 vr) d ( ? 2 + A 2) [( x _ ^2 + 7 2]l+! 


(A2) 


Further following Kong and Ravndal we set a = y 2 /(x — x 2 ) and introduce a second Feynman 
parameter to write 


h / 0 = 4yrail4 


2 e 


r(i + l) r 


dx- 


( 47 t) 2 T(2) Jo (x~x 2 )2 + 2 


x fto t ^ 


—- + 1 
U 2' 2 


0 J (2vr) d [( 1 _ a; )( g 2 + A 2) +a; ( g 2 + a )] 


3,6 
2 ' 2 


• (A3) 


Noting that 


(1 - c o)(q 2 + A 2 ) + uj(q 2 + a) = q 2 + (1 - ca)A 2 + u- 7 


x — x z 


(A4) 


we can carry out the integral and arrive at 

h/o = 4vraM 2 (^) 2C - r ^t^ r / 3 e) A f dx - fa 

n \ 2 J (4vr) 3 -r(2)r(| + f) J 0 ( x ~x 2 )^Jo 


u 2^2 


(1 — ta)A 2 + u 


7 


2 \« • 


x — x* 


(A5) 

Let us consider first the case A = 0, as done by Kong and Ravndal. In that case, the two 
integrals are 


dx- 


'0 [X — X 2 ) 2 + 2 Jo 


1 - 
U 2' 2 

dui 


= 7 2e I dx (x — x 2 ) 2 + 2 / dtata 2 2 


1 0 - 


7 


2 \ e 


x — X A 


= 7- 21 LNAEiliil _1_ ■ (A6) 


r(i + l) 1 


Hence, 


x= 0 _aM 2 N (n^Y [r(| + |)] 2 r(e) 1 


sir = 


8tt 3 /2 V 27 2 J r(2)r(| + f)r(i + |) 1 


47T 


i + 2 !og^-C E 

e 2 y 


+ 0(e). (A7) 
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Note that this differs slightly from the result obtained by Kong and Ravndal due to an error 5 
in Ref. [39]. If we keep the photon-mass term, the integrals over the Feynman parameters 
cannot be separated. After a couple of steps one arrives at 


dx- 


'0 ( X — X 2 ) 2 + 2 Jo 


dui 


U) 2^2 


(1 — u)\ 2 + LU 


= 2tt + 2n e 


7 


2 \ e 


X — X' 


1 + 2 log JJ- + 2 log (^1 + 

2y V 2 7 


+ 0(e 2 ), (A 8 ) 


so eventually we get 


aMl 

SIo = ^r 


+ 21 °s^- 2l °g( 1 + ^]- c '* 


+ OU) 


which reduces to Eq. (A7) when A = 0 


(A9) 


b. Simple momentum cutoff 


We now consider the integral with a simple momentum cutoff A. In that case, Eq. (Al) 
(with A = 0) becomes 


SlJ =0 (k) = 


d 3 qi f A d 3 q 2 


(27t) 3 J (27t) 3 k 2 — qf + ie ql k 2 — (qi + q 2 ) 2 + \£ ’ 
or, more conveniently for the present calculation , 6 


(A10) 


SlJ =0 (k) = AnaM 2 N 


d 3 gi f A d 3 q 2 


(27r) 3 ,/ (27 t) 3 k 2 — qf + is (qi — q 2) 2 k 2 — q 2 + ie ' 
This is logarithmically divergent, so based on dimensional analysis we know that 


(AH) 


SI A =°(k) = 


aM 2 N 
47r 


A log ^ + B + 0(k/A) 

rC 


(A12) 


with some dimensionless constants A and B. To determine these, we first carry out the 


angular integrals in Eq. (All) (all except one of which are trivial). As done in the calculation 
with dimensional regularization, we also analytically continue to k = by, 7 > 0 and let £ —> 0 . 
This gives 


sit°(k )=^ f <% q \ f d 92 21,192)1 . (Ais) 


7T 


'0 


' 0 


2gig 2 (g? + 7 2 )(?1 + 7 2 ) 


5 The argument of the Gamma function in the denominator of Eq. (100) in Ref. [JS] should be 3 — d/2 
instead of 2 — d/2. 

6 The shift involved in this rearrangement, which is equivalent to a different momentum assignment in the 
corresponding Feynman diagram, should only affect the terms suppressed by inverse powers of the cutoff. 



























Carrying out the integration over qi gives 


h/ 0 A =°(fc) 


ctMjj f A dqq 

2vr 2 J 0 q 2 + y 2 


log 


A ~q \ 
A + q) 


log 


A 2 + 7 2 \ 
q 2 + 7 2 / 


+ Li2 


g~A \ 

<1 - n) 


Iji-2 


q + A \ 

9 - 17 / 


+ Li2 


g~A \ 

9 + 17 / 


Li2 


, 

9 + 17/ 


(A14) 


where we have renamed the second loop momentum q 2 —> q. The remaining integral can be 
carried out using an indefinite integral which is given by a complicated sum of logarithms 
and poly logarithms. We do not write this out here, but simply note that the final result can 
be asymptotically expanded for large A to yield 


A = 1 and B = - ( log2 + , 


so that overall 


sit 0 = aM » 


° 47T 

where for convenience we have defined 


log -~C c + 0(k/ A) 
7 


Q- — log 2 + ~ 1.11943. 

Ati z 


(A15) 


(A16) 


(A17) 


Adding 0 = \og(aM N /2) — log(aM N /2 ) to Eq. (A16) we arrive at Eq. (41) 


2. Fully dressed bubble 


The fully dressed pp bubble shown in Fig. [3] resums the contribution of an arbitrary 
number of Coulomb-photon exchanges. We discuss it here in a number of schemes, starting 
with the results of Kong and Ravndal who used both dimensional regularization and a 
simple momentum cutoff. As discussed in the main text, Ref. [32] splits up the full bubble 
integral into a finite piece that gives the function H{rf) in the Coulomb-modified effective 
range expansion, and a divergent piece, which is the only integral that is regularized. In 
dimensional regularization, Kong and Ravndal find, for d = 3 — e, 


rdiv _ 

J 0 — 


Afjv ( ~ 


2vr d / 2 


(2n) d T(d/2) J 0 


d q q 


4-3 27r? /(?) 


e 2nri(q) _ 1 g2 


aM^ 

An 




where the result in the second line is obtained through the substitution x = 2nrj(q), which 
leads to the integral 


dx 


x 


e—1 


= r(e)C(e) • 


(A19) 


A subsequent expansion of the Zeta and Gamma functions gives the PDS pole in d — 2 
and the 1/e divergence for d —> 3, respectively. With a simple cutoff, the result is given in 
Eq. (1341) . 
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Note that Kong and Ravndal only consider the full series (empty bubble, one photon 
exchange, two photon exchanges, etc.) shown in Fig. [3] at once, although the logarithmic 
divergence really comes only from the single-photon contribution.' This leads to the puzzling 
result that while the 1/e divergences for the single-photon and fully dressed bubbles (see 
Eqs. (A7) and (A18), respectively) agree (same prefactor/residue), the logarithmic pieces 
that depend on the renormalization scale /i, as well as the constant terms, do not. Such 
an effect can be created by splitting up a formally divergent expression and applying the 
regularization only to a part of it. To see this issue, note for example that the integrand in 
Eq. (A18) cannot be expanded to order a to recover the single-photon result. 

this does not actually play a role- 


While in Ref. [39] this does not actually play a role—although they discuss the single¬ 
photon bubble for illustration, the authors eventually only use the fully dressed bubble for 
the matching to the Coulomb-modified scattering length—in the current work we require a 
scheme that is consistent throughout, i.e., we want to isolate the divergent single-photon 
piece within the fully dressed bubble. In Sec. A3 we first establish in configuration space that 
the logarithmically-divergent pieces should indeed be the same for both bubble diagrams, 
before we finally discuss in Sec. A 4 our approach of using the Coulomb T-matrix to isolate 
the the single-photon contribution, while still giving a closed expression for the rest. 


3. Configuration space 


Let us consider the two bubble integrals in configuration space in order to analyze their 
divergences more carefully. To this end, note that what we really want to calculate is the 
zero-to-zero Green’s functions [39] @0], 

Jo(k') = G C (E; r' = 0, r = 0) , E = k 2 /M N . (A20) 


Of course, this is infinite, so we want to isolate the divergence by starting from a finite 
separation and considering the case where it goes to zero. This procedure only deals with 
well-defined finite quantities, so we expect it to avoid subtleties associated with splitting up 
formally infinite terms. 

The Coulomb Green’s function satisfies the equation 


1 

2/i 



G c (B;r,r') = 6 <3 >(r - r'), 


(A21) 


in exact correspondence to the free case, 

^ (A r + k 2 ) G ( 0 +) (E- r, r') = 6® (r - r'). (A22) 

For the present case, the reduced mass is fi — M^/2. This Green’s function is known ana¬ 
lytically [72] [73] when one argument is set to zero. Adjusting to agree with the conventions 
used in Ref. m, we have 


G C (E- r' = 0,r) 


r(l - irj) W ir] .i/ 2 {-2ikr) 
47T r 



(A23) 


' In addition, the PDS pole in d, 


2 corresponds to the linear divergence of the empty bubble. 
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Expanding first the Whittaker function around r = 0, and subsequently the whole right-hand 
side in powers of a, we find 


G c (E-,r = 0,r) = 


M 


N 


47T 


—hi k + o/.Mn ^1 — Ce — log(—2iA;r) j 


+ Q{r ). (A24) 


The 1/r and log(r) terms in the above expression directly correspond to the linear and 
logarithmic divergences that we found previously in momentum space. 


The free Green’s function, ie., the solution to Eq. (A22), is given by 

-,ifc|r—r'| 


G‘ +) (i5;r,r') = (r|Gi +l (ii;)|r') = 


(+), 


/i e 


27T |r — r' 


(A25) 


The first two terms in Eq. (A24) come from, with // = Mjv/2, 

M n 


G< +, (S;r' = 0.r) = 


47T 


1 .. 
—hi k 
r 


+ 0(r ) 


(A26) 


This corresponds to the contribution of the empty bubble. 

We now turn to the single-photon diagram and expect to find exactly the remaining 
divergence, because diagrams with two and more photons should be finite. What we have 
to calculate is thus 

G ( c\E-y, r) = {r|G< +) (B)CcGh , (S)|r'). (A27) 

where 

(rlVd'lr') = (r — r') — (A 28 ) 

r 

is the Coulomb potential. Setting k = + \e and r' = 0 in Eq. (A27) gives 


G$(E- 0,r) = 


~2V 


2 roo 


dr" e ikr " 


dcos 6 




3 ifc|r—r"| 




, r" ■ r = r" r cos 6 , 


(A29) 


where we have inserted identities to separate the operators and carried out the trivial in¬ 
tegrals. The angular integral simply recovers the well-known expression for the S- wave 
projected free Green’s function: 


dcos# 


e ifc |r-r"l 2 e ikr > sin (kr<) 


i -i 




k 


, r > — max(r, r") , r< = min(r, r"). (A30) 


The remaining integral in Eq. (A29) can then be done by splitting up the domain according 
to 


dr" = / dr" + 


dr". 


(A31) 


In each region, r< and r> are then uniquely defined. The first integral is 

sin(fcr /; 


rr gihr" 


-dr" = ^ [c E + log(2fcr) — Ci(2fcr) — i Si(2 kr)j , (A32) 
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in terms of the sine and cosine integrals m, respectively 


Si(-) 

Ci(s) 


sin(f) 


d t = 0(z ), 


cos(f) 


dt = C E + log z + 0(z 2 ). 


(A33a) 

(A33b) 


For the second integral we find 



J2\kr" 


-dr' 1 


r(—2i kr). 


(A34) 


Combining the results, we get 


G') ) (£;r' = 0,r) = 


a/j, I le 


7T kr 


-Akr r 


Ce + log(2 kr) — Ci(2 kr) — i Si(2fcr) + sin(fcr)r(— 2ikr) 

(A35) 


Finally, expanding this around r = 0 and inserting /j = Mat/ 2, we arrive at 

1 — C E — log(—2i kr) + 0{r). 


G^\E- r' = 0,r) = 


aM 2 N 

47r 


(A36) 


This is exactly the term linear in a in Eq. (A24). 


4. Fully dressed bubble with T-matrix expansion 


Having established that the remaining divergence indeed comes exactly from the single¬ 
photon diagram, let us take another shot at consistently isolating it in a momentum-space 
calculation. To this end, note that the Coulomb Green’s function satisfies the operator 
equation 

G C (E) = G ( 0 + \E ) + G ( 0 + \E) V c Gc(E) . (A37) 

On the other hand, we have the Coulomb T-matrix defined by the Lippmann-Schwinger 
equation 

T C (E) = V C + V C G ( o + \E ) T C {E). (A38) 

With 

VcGc(E) = T C {E) Gq + \e) , (A39) 


we can write Eq. (A37) as 

G C (E) = G< +) (£) + G‘ + >(£) T c G { +\e) . 


(A40) 


The advantage of this operation becomes apparent once we note that the Coulomb T-matrix 
is known in closed form. For example, in terms of the Coulomb potential in momentum space 


Vc(p,q) = 


47 ra 


p -q) 2 ’ 


(A41) 
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one has 


T C (k; p, q) = (p|Tc|q) = V c (p, q) < 1 - 2 ip 


s + 1 
s — 1 


- 1 V 


ds 


1 — e 


(A42) 


or 


T c (/c; p, q) = V c { P, q) < 1 - A 2 F\ 1, ip, 1 + ip; 


A - 1 
A + 1 


- 2-Pi 1, if 7, 1 + i? /; 


A + 1 
A - 1 


, (A43) 


with 


e = 


(p 2 — kr){q 2 — k 2 ) 


A 2 = 1 + e. 


(A44) 


fc 2 (p — q) 2 

For these and a number of alternative representations, see for example Ref. m Note that 
both Eq. (A42) and (A43) express the T-matrix as the Coulomb potential plus a remainder 
that summarizes the exchange of two and more photons, 



T C (k ; p, q) = V c (p, q) + T£(k] q, p) . 

(A45) 

Hence, we 

can write the fully resummed bubble integral as 



Jo(k) = I 0 (k) + 5I 0 (k) + 6J 0 (k), 

(A46) 

where 

h(k) = M N f f q - 1 

J (27t) 3 k 2 — q 2, + \e 

(A47) 


is the empty bubble, SI 0 (k) is exactly the single-photon bubble calculated in Sec. Al, and 
5Jo(k) remains to be investigated. The crucial point is that it is supposed to be ultraviolet 
finite. 

The integral we have to evaluate is 


SJ 0 (k) 


r2 [ d 3 P [ d 3 g _ T c(k\ g,p) _ 

N J (27t) 3 J (27t) 3 (p 2 — k 2 — i e)(q 2 — k 2 — ie) 


(A48) 


Noting that T^{k] q, p) depends only on the angle between q and p, we can carry out 
all but one angular integral. The remaining one over cos 9 is then an S'-wave projection. 
Fortunately, partial-wave projections of the full Coulomb T-matrix are known in closed 
form jZSl mi- The S'-wave result can be written as m 


Tcp(E = k 2 /(2fj.)\p,q) 


where 


in k \ . 

-2-F1 (1,177 ,1 +W,aa) - 2 

p pq L 

+ 2^(1, i/7,1 + i?7; l/(a'a)) 


Fi(l, ip, 1 + i 77 ; a/a') 

- 2 -Fi(l,i? 7 , l + ipja'/a) , 


p — k , q — k 

a = - - , a = - - . 

p + k q + k 


(A49) 

(A50) 
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With the projection of the Coulomb potential alone, 


Vc,o (p,q) = —Qo 


pq 


p 1 2 3 4 5 6 7 + q 2 
2 pq 


it a 

= —log 
pq 


2 pq + p 2 + q 2 
2 pq — p 2 — q- 


(A51) 


we can write 


SJo(k ) = 


M 


N 


Y 


27T 2 J 


d pp 2 dqq 


2 T c ,o(E = k 2 /(2p);p,q ) - V c ,o(p,q) 
(p 2 — k 2 — E)(g 2 — k 2 — is) 


(A52) 


This integral looks difficult to solve, but at the very least it can be evaluated numerically. 
We do this with a (large) momentum cutoff and a small imaginary part of k in place for 
regularization. Considering the real part, one finds that it is consistent with 


Re5J 0 (k) = 


aM 2 


N 


47r 


(. 


Re jV’(ih) + 2 i^ f + 


with a constant 


C A « 0.579. 


(A53) 


(A54) 


The level of agreement is demonstrated in Fig. [l4](a). Note that Eq. (A53) contains the 
function H(p) except for the log ir/ part. This agrees with what one would expect here 
because in the sum (A46) we already have a log k term from the single-photon contribution 
Slo(k). Considering the imaginary part, we find that it agrees very well with 


ImhJo(fc) = 


aM 


N 


47T 


Im 




M^k 
47r 


(A55) 


as shown in Fig. 14(b). Altogether, we have thus established numerically that 

SJo{k ) = 


aM 2 N 
47r 


^ + l, + c *) + A t ik - 


(A56) 


We expect that increased numerical accuracy would yield Ca —>• Ce ~ 0.577216. The value 
in Eq. (A54) represents our most accurate result (obtained with a momentum cutoff of 
16000 MeV and Irnfc = 0.05 MeV); decreasing the former and increasing the latter gives 
larger values for C'a, so indeed we see a convergence pattern. 


[1] H. A. Bethe, Phys. Rev. 76 (1949) 38. 

[2] H. A. Bethe and R. Peierls, Proc. Roy. Soc. A148 (1935) 146; ibid. A149 (1935) 176. 

[3] P. F. Bedaque and U. van Kolck, Phys. Lett. B 428 (1998) 221 [nucl-th/9710073]. 

[4] U. van Kolck, Lect. Notes Phys. 513 (1998) 62 [hep-ph/9711222]. 

[5] D. B. Kaplan, M. J. Savage and M. B. Wise, Phys. Lett. B 424 (1998) 390 [nucl-th/9801034], 

[6] P. F. Bedaque, H.-W. Hammer and U. van Kolck, Phys. Rev. C 58 (1998) 641 [nucl- 
th/9802057], 

[7] D. B. Kaplan, M. J. Savage and M. B. Wise, Nucl. Phys. B 534, (1998) 329 [arXivmucl- 
th/9802075]. 


34 


















Figure 14: Numerical evaluation of 8 Jq( k): (a) real part; (b) imaginary part. Black dots correspond 
to the integral given in Eq. (A52). The solid lines show the functions defined in Eqs. (A53) (a) 
and (A55) (b). The calculation was performed with an explicit small imaginary part of e = 
0.05 MeV added to the real momentum k and a momentum cutoff of 16000 MeV. 
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